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Abstract 

We investigate the effect of a constant external non-Abelian field on chiral 
symmetry breaking in a (2+l)-dimensional Nambu-Jona-Lasinio model and 
in 3D QCD by solving the gap equation and the Bethe-Salpeter equation, and 
also by RG analysis. In the (2+l)-dimensional NJL model chiral symmetry 
breaking occurs for any weak coupling constant but in 3D QCD catalysis of 
chiral symmetry breaking does not occur. 
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Though there are increasing evidences for chiral symmetry breaking (%SB) in QCD [|TJ, we 
still do not understand clearly how it occurs, due to a lack of comprehensive understanding 
the QCD vacuum. One of the popular analytic methods to study %SB in QCD is by solving 
the Schwinger-Dyson (SD) equation for the Green functions, from whose solution one can 
extract the structure of the vacuum 0. Most attention in SD analysis has been focused 
so far on the effect of the strong coupling in QCD; if the coupling is bigger than a critical 
coupling, the SD equation admits stable xSB solutions. But, the effect of self- interacting 
gluons in SD analysis has not been seriously considered except in the running coupling, 
even though it has been conjectured that the monopole condensate in QCD leads to color 
confinement and thus %SB 0, which is now shown rigorously in N = 2 supersymmetric 
QCD [§]. 

Since the quark condensate under an external gluon field A ^ is related to the density of 
eigenvalue of Dirac operator ip(A) at zero eigenvalue ||, 

(0|gg|0) A = -vr Pyl (0), (1) 

it is conceivable that a certain gluon field configuration which gives a constant eigenvalue 
density at small eigenvalue like a free (0+l)-dimensional spectrum is a source of %SB. On 
the other hand, it has been recently found that an external magnetic field catalyzes chiral 
symmetry breaking in the Nambu-Jona-Lasinio (NJL) model || and in quantum electrody- 
namics • The fact that a constant magnetic field enhances fermion dynamical masses 
was known from the papers fL0| . Studying (2+l)-dimensional NJL model with discrete 
chiral symmetry, it was shown |11|] also that chiral symmetry turns out spontaneously bro- 
ken for any B ^ 0. However, what was shown in || is not just that the magnetic field 
enhances dynamical fermion masses created by supercritical NJL interaction, but that it 
catalyzes generating the mass even at the weakest attractive interaction and this effect is 
universal, i.e., model independent in 3+1 and 2+1 dimensions^. The physical reason for the 
magnetic catalysis of xSB is that, in the presence of external magnetic field B = Bz, dimen- 
sional reduction occurs in fermion motion. The fermion spectrum becomes one- dimensional, 



E = ±|p z |, at energy much lower than the Landau gap E <C J \eB\ and, for one-dimensional 



fermions, an arbitrarily weak attraction at scale /i = \J\eB\ leads to a strong attraction in 
the infrared region ||. It is important that in |||7| continuous chiral symmetry was consid- 
ered, that is far nontrivial problem since the appearance of Nambu-Goldstone bosons along 
with reduction toD = + lorD = 1 + 1 dimensions might wash out symmetry breaking, as 



happens, for example, in the chiral (l+l)-dimensional Gross- Neveu model [13|. However, it 
was shown that notwithstanding dimensional reduction Goldstone boson propagators keep 
their three or four dimensional form. 

In studying %SB in QCD, it is natural to look for a gluon field that mimics a magnetic 
field and gives a (l+l)-dimensional spectrum for the Dirac operator. We consider this in a 
(2+l)-dimensional NJL model and in a (2+l)-dimensional SU(2) gauge theory for simplicity. 
These theories in 2+1 dimensions have been studied intensively recently, since they can serve 
not only as theoretical laboratories for investigating aspects of chiral symmetry breaking in 



§For recent application of this effect in condensed matter physics see pi 
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QCD but also describe the high temperature thermodynamics of 4-dimensional theories as 
well as certain planar condensed matter systems. 

A constant and uniform chromomagnetic field for a non-Abelian gauge theory can be 
described by two different types of gauge potentials jl4yi5|| , a feature not available for the 
abelian case. As is known these two different potentials lead to two different energy spectra 
for fermions in such a field [14], 16]. The non-Abelian external field which mimics a mag- 
netic field and has the largest global symmetry including explicit translational invariance, 
is described by a constant gauge potential 



.4, 



2 



with A\ = A\ = JH/g, and others = 0, 



(2) 



where r a 's are the Pauli matrices (a = 1, 2, 3), g is the gauge coupling constant, and if is a 
constant. Then, the field strength tensor, 



(3) 



so that F? 



i 2 — — F21 = —H and others are zero, where the covariant derivative 



dp + igA a T a /2. Under the external field, Eq. @, the fermion propagator satisfies 

\f {id^ - gA^) - m] S(x -y) =i6(x-y), (4) 

where m is a fermion mass and we take the four- component reducible (2+l)-dimensional 
spinor for fermions. To find the fermion propagator we look for a solution of the form 



S(x -y) = i [7 M (ify - gAp) + m] A(x - y), 
where the Fourier transform of A (re — y) satisfies the following equation: 



P P^ 



m 



A(p) 



where 

With our choice of the vector potential 



a 



{Y,Y} = 2g^. 
the equation (pi) simplifies to 



P 



2h(p- f) - 2K 1 + 2/iW 2 



m 



A(p) = 1, 



(5) 



(6) 



(7) 



(8) 



where we introduced the notation h 
we find: 



1 gH /2. The last equation can easily be solved and 



A(p) 



1 



Q1Q2 



P 



2^3 12 



2h z + 2hp- t- 2h z r 6 a 



m 



where 



QiM =pl~m 2 - 



p 2 + h 2 ± h 



(9) 



(10) 
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Finally, the Fourier transform of the fermion propagator reads 



S(p) = i 7 M P M + m A(p). 
We see that if m = the pole of the propagator occurs at 



(11) 



p = ±(J I P + h^±h). 



(12) 



At low momenta, \p\ <C h, the spectrum changes drastically, it has a non-relativistic branch: 
Po = ±pP/(2h). It occurs that just this mode is responsible for chiral symmetry breaking in 
NJL model. 

Now, we consider the (2+l)-dimensional NJL model with N four- component fermions in 
the presence of the above non-Abelian field, described by the Lagrangian density 



2 + f^75^) 2 + (i>l3^ 2 



(13) 



which has U(2N) symmetry. In terms of auxiliary fields a, r and 7r it can be rewritten as 



C' = ipi^D^ip - ip (a + 7 3 r + rfir) ip - ^- (o 2 + vr 2 + r 2 ) . 



2G 



(14) 



In the large N limit the vacuum is determined by the stationary point of the effective action, 
obtained by integrating over the fermions: 



JS cB (ct,t,tt) 



i J d 3 xC 



(15) 



or 



N r 

SeS = ~2GJ + n2 + T ^ ~ %Tr hg ^ D ^ ~ ( a + i757r + 73T ^- 



(16) 



Since the vacuum should respect translational invariance, we need to calculate the effective 
action for constant auxiliary fields to find the vacuum. In this case, the effective action is 
just S e s(a, t, 7r) = — V(a, r, n)L 2 T, where L 2 T is the spacetime volume and V is the effective 
potential. The potential V depends only on the chiral invariant p 2 = a 2 + n 2 + r 2 and for 
its calculation it is sufficient to consider a configuration with r = 7r = and a = const. In 
the proper-time formalism we get for V(cr): 



V(a) 



Na 2 



f ds f d 3 p isl p u pn_3 ■ (T »u F « vT a_ a 2} 

J s J (2n)3 



2G 2 J s J (2tt 



(17) 



The matrix in the exponent of ( |17D has two eigenvalues (compare with (|T0|)) 



\ 1 ,2=p 2 o-v 2 -WP 2 + h*±h} 



(18) 



each of which has 4N degeneracy. Thus we write 
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V(a) 



Na 2 
~2G 



+ 2N 



d 3 p 



l/A 2 



ds 

~7 J (2tt) 3 



-s\p 2 3 +a 2 + (^p2+h^h) 2 } 



+ e 



s\pl+a 2 + {^P 2 +h 2 +h?\ _ 2e -M+P 2 ) 



(19) 



where we made a rotation to Euclidean space, subtracted the part corresponding to free 
fermions (a = 0, h = 0) and introduced the ultraviolet cutoff A. The effective potential up 
to the order of A is computed to be 



V(p) 



Np 2 N 
~2G + 7 



p 2 A , p 3 + (p 2 + Ah 2 f 2 u2 nr — 2 

+ - ^ h 2 Jp 2 + Ah 2 - In 

"K 6 v I 



p 2 h^ (2h + ^p 2 + Ah r 



(20) 



where we substituted a by the chiral invariant p. The expression ( p0|) agrees with that one 



derived in |T1J, except the SU(3) color group was considered there. The stationary condition 
for the effective potential gives the gap equation: 



Np N 



2pA 



p 2 + p (p 2 + Ah 2 ) 1 ' 2 - ph In 



2fe + Vp 2 + Ah 2 
P 



As /i —>■ 0, we recover the known gap equation 

p 2 = pA / ' 



7T 



— V 

2AG 



(21) 



(22) 



It admits a nontrivial solution only if the coupling G is supercritical, G > G c 



TV 



3/2 



/2A. 



The magnetic field changes the situation: at h ^ 0, it is easy to see that a nontrivial solution 
exists at all G > 0. The reason for this is that the magnetic field enhances the interaction 
in the infrared region: at h ^ the last term in fl2lD is responsible for the existence of a 
solution at any G > 0. The nontrivial solution always minimizes the effective potential, if 
the four- Fermi coupling G > 0, and therefore is the vacuum solution. Thus, unlike the NJL 
model with h = 0, chiral symmetry breaking occurs for any value of G > 0. In particular, 
when G — > 0, 



p = Ah exp 



7T 

hG 



(23) 



while p ~ O.Ah when G ^ G c , the critical coupling for chiral symmetry breaking of the 
NJL model in the absence of the external field. Since p defines the fermion mass in the 
Lagrangian density (0), we find that the fermions get a dynamical mass, which is m^ jn = 
Ah exp [— 7r/ (hG)] for weak coupling. To analyze the supercritical region G ^ G c , it is useful 
to introduce the dynamical mass m , defined by the Eq.(^). For weak fields (h <C m ), we 
find 



"idyn = p = m 1 + 



3h 2 
AttIq 



(24) 



i.e. mdyn increases with h. 
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If the cutoff A = h, the calculations get much simpler and the reason for the catalysis 
is easy to understand. Let us now introduce a dimensionless coupling A = hG. Then the 
effective potential becomes (p <C h) 



V(p) = Nh 



2h 2 p 2 1 „/ p 2 4 
3tt 2A 4vr \ 16h 2 ^ 



7T 



0(1), 



(25) 



from which we can read off the scaling behavior of the coupling G by taking a different 
renormalization point as 



(Nh) 



.1 d 2 V 



dp 1 



1,1, jfL 

A 2?r I n 16/i 2 



7T 



+ 6 



(26) 



The dynamical mass of a fermion, which will be invariant under the change of the renormal- 
ization point, is of the order of 



/ 7T 2 

4/iexp — — — + 7r H — — - 3 
V hG \Pk , 



(27) 



This result is very similar to magnetic catalysis, where an arbitrarily weak magnetic field 
catalyzes chiral symmetry breaking, though the motion of fermions at low energy is not 
dynamically reduced to one dimension by the uniform chromomagnetic field. The reason 
for catalysis of by the uniform chromomagnetic field is that the spectrum of fermions 
changes at low energies, E <C h, so that the scaling dimension of fermions is the same as 
that of fermions in (l+l)-dimensions and thus the four-Fermi interaction becomes a relevant 
operator. This can be easily seen by performing the renormalization group analysis similar 
to that which was applied to the ordinary magnetic catalysis M. 
Let us expand the fermion field as 



W,t) = / J2Mt,k)e ikr u A {k), (28) 
Jk A 

where UA0) are the eigenfunctions of the operator a ■ (k — gA a; Y^j with eigenvalues Ea = 



(i i v k 2 + h 2 + (3hj each of which is doubly degenerate. (A denotes collectively a, (3 — ± 
After Fourier-transforming in t, we find that the kinetic term in the action is 

du d 2 k ~ 



s = E 



2tt (2tt) 



ip\i(u, k)(u - E A )ilJAi{u) 



1,2, 



(29) 



where ipAii^, k) is the Fourier transform of ipAi{t, k) and the spinors um are chosen to be 
normalized as u Ai UAj = &AB&ij- 

When \k\ h, the effect of the external field is important. It changes the energy 
spectrum drastically. To derive the low energy effective action, we integrate out the modes 
with \ k\ > h. Then, the kinetic term in the effective action becomes 



s = E 



du d 2 k 
\U\<h 2^(2vr) 2 



k 



+ 



(30) 
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where we have expanded the energy eigenvalues Ea in powers of momentum and the elipsis 
denotes the terms with higher powers of momenta. 

To determine the scaling dimensions of the fermion modes, consider a transformation, 



UJ — > SLd, 



k -> s 1/2 k with s < 1. (31) 



Since the kinetic term must be invariant under the scaling transformation, we find ip a i- has 
scaling dimension —3/2, while ip a i,+ has —1. Namely, in the coordinate space, V'm.iC^j t), the 
Fourier transform of ip a i,±, has scaling dimension 1/2 and 1, respectively. We see that the 
component ipad- is more relevant in the infrared region than the component if) a i,+- Because of 
the change in the energy spectrum due to the external field, the scaling dimension of fermion 
field changes for the component ipai- at low energy as if the spacetime is dimensionally 
reduced from 2 + 1 to 1 + 1. Because of this dimensional reduction, the four- Fermi interaction 
of ifjai- is marginal. Now, let us further integrate out the modes of sh < uj < h, s^ 2 h < 
\k\ < h to find the change in the four- Fermi interaction. The one- loop correction to the 
four- Fermi coupling in the leading order in 1/N is given then 

5G = ^-[ Tr[S(k)S(k)] = -—hi8. (32) 

iV J S h<U<h,S l / 2 h<\k\<h 71 

In terms of the dimensionless coupling A, we find 

= (33) 

Therefore, the four-Fermi interaction becomes a relevant operator at quantum level if G > 0, 
which at low energy leads to strong attraction between fermion and antifermion to form a 
condensate. The dynamical fermion mass, the RG invariant scale, is m^ yn — hexp[—ir/hG)], 
which is, in leading order in the coupling expansion, the same as the dynamical mass we 
obtained by solving the gap equation fl2"T|) for weak coupling A at \i — h, where the RG 
analysis is reliable. 

The (2+l)-dimensional SU(2) gauge theory is more complicated than the 3D NJL model, 
since not only the fermions but also the gauge field propagate nontrivially under the external 
uniform chromomagnetic field. The gauge field propagator (in Feynman gauge) in momen- 
tum space is determined as a solution to the equation: 

VpVgap - 2zgF^I a ] S»(p) = (34) 

where = — gI a A a ^, and I a are the generators in the adjoint representation. As is seen, 
to solve this equation one has to invert a matrix with double indices, Lorentz and internal. 
The structure in Lorentz indices, however, can be simplified considerably. That is due to 
the fact that the only nondiagonal Lorentz matrix appearing in the equation is F^g, and it 
commutes with g a/3 . So, we represent the gluon propagator through the projectors which 
diagonalize the F^-matrix: 

S fh(p) = -i £ S j (p)A^ , (35) 
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where 



p3(3vp3 7 



H 2 



£ A 



Pi 

3 



■37 



A 



±i 



H 



F 3f3 Av-y 



iRjAf. 



H 2 



(36) 
(37) 



Thus, we obtain the equation for 



V 



Ah 2 ((I 1 ) 2 + (I 2 ) 2 ) - 4h(p- f) + 8jh 2 I 3 



S 3 = 1, 



which is easily solved by a simple inversion of 3 x 3 matrices. The result reads: 



S J 



Uj + b(pl) + c{Pf + jdP + f (pi) - 2hjl c 



E; 



where 



aj 



(p 2 - Ah 2 ) 2 - 16/iV - Q4h 4 j 2 , b = 4h(p 2 - Ah 2 ), 
-4h 2 (p 2 -4h 2 ), d= -8h 2 (p 2 -8h 2 ), f = lQh 2 , 
- [p 2 - 8h 2 ){p 2 - Ah 2 ) 2 - 16hY{p 2 - Ah 2 ) - 6Ah*f{p 2 - 8h 2 ) 



(38) 



(39) 



(40) 
(41) 
(42) 



We find that under the external chromomagnetic field the gauge field is not massless, 
instead its propagator has poles at nonzero p, all of order h, satisfying (p 2 — 8h 2 )(p 2 — Ah 2 ) 2 — 
IQtffPij) 2 — Ah 2 ) — 6Ah 4 j 2 (p 2 — 8h 2 ) = 0, j = 0, ±1. Seven of these poles correspond to 
massive modes while the remaining two describe tachyonic modes signaling instability of the 
perturbative vacuum in the external chromomagnetic field (g) (analog of the Nielsen-Olesen 
instability in four- dimensional Yang-Mills theory for abelian configurations of chromomag- 
netic field ||18||).p*l Since we are interested in infrared dynamics of fermions with energy and 
momenta E, \p\ -C h, the presence of tachyonic modes of order h is not very important. In 
fact we will be working in an effective low energy theory where h serves as ultraviolet cutoff. 
Further we need the gauge propagator which is approximately constant for E,\p\ ~ 0, 



3a r3b 



6h 2 H 



Eaf3 ^3ab _|_ 



(F 2 



la/3 



3h 2 H 2 



(S<* - 5 3a 5 3b ) . 



(43) 



Since the gauge field is massive, the gauge field exchange interaction will generate a four- 
Fermi interaction at low energy E <ti h: 



4 {x) 



A8h 2 



+ ■ 



(44) 



where the ellipsis denotes other than NJL interaction terms. For a while we assume that only 
NJL-like interactions play important role in the generation of a dynamical mass. Then, the 



**For a thorough investigation of the instability of non- Abelian fields given by constant potentials 
see @. 
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„2 



effective four- Fermi interaction generated by the massive gauge field is repulsive (— ^2 < 0). 
Therefore, even though the four-Fermi interaction for the ip a i t - component of the fermion 
is marginal as we have shown in RG analysis of the 3D NJL model, its /3-function is pos- 
itive and it becomes an irrelevant operator at one-loop. Thus we may conclude that the 
chromomagnetic field does not induce chiral symmetry breaking in 3D QCD. 

Since, while deriving the effective four-fermion interaction, we omitted terms other than 
those in ([44) in order to apply the RG analysis, below we confirm our conclusion by analysis 



of the gap equation for a fermion dynamical mass. There are two equivalent approaches 
one can follow in studying mass generation, that of Schwinger-Dyson equations ||||, and 
that of Bethe-Salpeter (BS) equation M. We consider here the BS equation for a gapless 
Nambu-Goldstone (NG) boson, composed of a fermion and antifermion, which must appear 
when the chiral symmetry is spontaneously broken. 

The homogeneous BS equation for the NG bound state reads: 

P = "to / ^jrS(p)pS(p)^S;l, (45) 

where a = (? 2 / (47r), and the Bethe-Salpeter wave function for the bound state is given in 
general as 

p = 7 5 (A + r 3 a 12 B + (7 ■ f)C) . (46) 

Since the dynamics responsible for fermion mass generation is due to momenta near p = 0, 
we can neglect the momentum dependence in the gauge field propagator in the BS equation. 
Thus, the functions A, B, C become momentum independent and Eq.(|45|) turns out to be 
an algebraic system of equations. 

By direct substitution of (^6|) in the BS equation, we observe that the equation for 
coefficient C decouples from the set of equations for A and B. The latter is: 

A _ma r d 3 p [6h 2 + 13(p 2 - m 2 )] A+ [6h 2 - 16(p 2 + 2f -m 2 )]B 
ma r d 3 p [78h 2 - lQ(p 2 - m 2 )} A + [78h 2 - 23(p 2 + 2f - m 2 )} B 

~2^J WMp) ' (48) 

In the limit m — > 0, this reduces to the following 

B^~^(^)(A + B), (50) 



or in other words, we obtain 



7 "ln(A). (51) 



32h \ m 



Thus, there is no solution satisfying the condition m <S h, which agrees with the result of 
RG analysis (analysis of an equation for C-coefficient shows no nontrivial solution, too). 
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While 3D QCD itself is known to be unstable in perturbation theory ||19 , and might lead 

the presence of an exter- 



to generating a mass for Yang-Mills fields nonperturbatively |20 



nal chromomagnetic field enhances this instability. A chromomagnetic field influences the 
behavior of charged vector particles more strongly than that of fermions and further study 
is needed to make decisive conclusion about the magnetic catalysis in 3D QCD. Perhaps, si- 
multaneous study of generating chiral condensate and a condensate of vector particles would 
shed light on this problem. The absence of chiral condensate in a chromomagnetic field in 
3D QCD may cast doubts upon straightforward application of the results of a magnetic 
catalysis in QED4 M to QCD4 El. 



In conclusion, we investigated the effect of a constant external non-Abelian field on chiral 
symmetry breaking in a (2+1) -dimensional Nambu-Jona-Lasinio model and in 3D QCD. We 
studied the effect by gap equation analysis and RG analysis and found that in the (2+1)- 
dimensional NJL model chiral symmetry breaking occurs for any weak coupling constant 
(in accordance with earlier findings JHJ) but catalysis of chiral symmetry breaking does not 
occur in 3D QCD. 
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